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Introduction 


In connection with Project FLAMBEAU it was thought useful to make a 
collection of simple theories relating to mass fires. The list is far 
fran exhaustive; in particular attention has been concentrated on the 
conditions at peak, with ignition, growth and spread being neglected. 
Plume theories have also been neglected, as being too complex; 

Thomas' Working Notes on Project Flambeau (1) give a survey of them. 

Null hypotheses have been stressed. The prediction of burning rate in a 
Flambeau type situation is regarded as the central problem, although later 
sections do cover theories related to shelter and measurements. 


The Descriptive Model-Clive Countryman (2) 
1 In temns of height: 


Fuel Zone 

Combustion zone-up to perhaps 100! 

Transition zone - with sane down drafts - up to at most 2! - 
this is the upper limit of downward influence 

Thermal convection zone - up to 1000! - 20,000! 

Smoke fallout zone - a thin layer 

Condensation convection zone. 
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Transition Zone 


Combustion Zone 


A Nig iy a Fuel Zone 


Theories of pile interaction 
2 At peak intensity, piles do not interact. 


2a The burning rate of a pile (b) is not affected by other piles, and is 
proportional to the surface area of fuel exposed. 


ALL BURNING AT THE SAME RATE 


[ 
2b An isolated pile exposed to standard met. conditions will burn at a 
rate b (0,t) varying with time t. The reasons for this variation are the 


initial spread of fire through the pile and the variation in the surface 
area exposed as the smaller twigs are consumed. 


FIRE SPREADING FINES BURNING OUT 


b is better expressed as a function of the total amount burmmt:~- 
t 
I =f b dt. 
[e) 


Thus the standard burning rate becames b (0,1). If the actual burning 
rate is b(t), then part of its variation is due to the above causes. 
Put b'(t) = b(t)/b (0,1) 


Then the variation of b'(t) with time is mainly due to variation of 
interaction effects. 


2c Pile interaction is by ignition or re-ignition (i.e. spread from 
pile to pile), in the case of close spacing under poor burning conditions. 


2d According to Thamas (3) the burning rate of a crib or pile in still 
air is proportional to 20-6 (A, yee , 


where A, = Surface area of fuel A, = Total vertical vent area 
Z = Height of pile 


A CRIB 


Ag: Surface area 


of fuel a Ay = Total vertic 


vent area 
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5) Piles produce their effects on the environment independently 


3a The environment at ea point is the sum of that. from each of the piles, 
each*»t the appropriate distance and burning st the appropriate rete. 


PLAN OF PILES 


3b The contribution from a pile is proportional to its burning rate. 


3c The contribution from a pile at distance d is proportional to 1/a? 


4 Pile interaction is by means of the wind generated. 


4a A single pile burning at a rate b produces a wind b.w(d) directed 
towards itself at a point distance d. In an array a pile has acting on 
it a total wind W «= Wo+ sum for other piles (b.w) 


Where Wo is the ambient wind. 


=w 
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PLAN OF PILES Ne baw(d2) 


Now the burning rate b depends on W :- 
be b(IWI) 
Hence feedback 


Moreover, if at a pile directly downwind of the array W is acting inwards 
this is a non-spreading fire; if outwards a spreading one. 


4b As theory 4a, but only the nearest ring of piles affect sum (b.w); 


all the rest are shielded. Hence only piles on the edge of the array burn 
at a rate different from b(Wo). 


iw 


PLAN OF PILES 


No Effect On X 


by b, byw(d)) x 


4e As theory 4b, but W 
at a rate b(0). = 


Wo 


No effect on X 
PLAN OF PILES 


Noeffect on X Dw 
zy bw(d) = 
ee 


4d As theory 4a, but determine w(d) experimentally in the presence of 
inert obstacles representing the other piles (in practice the actual other 
piles might be used, but not ignited and kept well damped down). 


4e As in theories 4a - 4d, but allow also for the variance of the wind and 


its effect on the burning rate. Assume that the variance of a wind can- 
ponent b.w is u b Iwl , where u is a constant 


10) 
Also b= b (IWI ) can be written as:- 


DEST DN(We) 


The total variance = Ge u ( + sum for other piles (b.iwi ») 


Now expectation of Wee wW +G 


Root Mean Square W 


This determines b 


TIME ——> 


so put b = b (W + G). 
4f The total angular mamentw is zero; if there is, say, a general 


clockwise swirl about the array, its angular momentum will be balanced by 
anticlockwise fire whirls 


Clackwise ratatia 


of main column 
Anticlockwise rotation 


of fire-whirl 


by When a proportion of the draught for one pile is cambustion products 
from the others, this will affect the fire because of the decreased oxygen 
and the increased temperature 
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Altered burning rate Standard burning rate 


5a Low level air up to some height h is readily available in the central 
area of a mass fire, If this is sufficient to burn up all the fyel there 
will be no important interaction. With a fire loading 1 cals/am™ needing 
a column of air 1 cms. high, this means 1 « h; if not vertical convective 


circulation will set in bringing air from above. 


ee “A= Saas 55 oe 


No interaction until shaded air has been exhausted 


5b With an ambient wind spoed No and a width of array L the volume of air 
available/sec is Wo Lh. 


Hence the burning rate is khLW. cals/sec. and thus the rate of advance 
of the fire front will be khLW)/1L &ns/sec = kh) /1 cms/sec. 


ca 


5c¢ As in theory 5b, let air concentration (fraction of fresh air by 
volume) be c, the burning rate/unit area R(c) cals/am2/sec., and the 
burning efficiency e(c). Then 1 c.c. of air, given the fuel, will give 
ke(c) cals. (2(1) = 1). Consider the area near the upwind edge of the 
array, where the ambient wind is the dominant air supply, and assume 
equilibrium conditions, so that R does not vary with time. Then heat 
generated/unit area = R = hw, (de) e(c) k where x = distance from upwind 
edge of array. dx 


therefore dc = -R(c) 
ax hWke(c ) 


Solve this for c(x) and hence find R(x) and e(x) 


In particular, for x = 0, (de) = oO) 
(dx) 9 hWk 


t21,28&3 


Sd Assume air is only available by vertical convection (see theory 5a). 


If it is available up to height h,, then for a vertical prism of 
air of height h, and unit cross-sectional area, heat generated/sec = 


h,(-de) k e(c) & R anddc _ site) 
dt dt hy ke c) 


As before, this can be solved for c(t) and hence for R(t) and e(t). 


R(x) 


5e As in theory 5d, but assume h 


1 is proportional to the burning rate, 
or h, = hk (hy is a constant). 


b Doubled 


Hence as R builds up to a maximm, more and more fresh air is 
entrained at the top, While this is happening, 


ch, + 1.dh, = (c + de) (hy rs dh, ) 


therefore dh, (1 - ce) 


(als ee 
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and an extra term must be added to dc, giving 
dt 


(1 - ¢) 


ac 1 
dt 


After the maximum there is no further entrainment and the extra term 
should not be added, 


hence dc __-1 
dt hoke(c) 


5f As above, but h, is proportional to the flame height F. 


As above with h, 


| proportional to 
} hn 


flame height 


b SMALL b DOUBLED 


If theory 7b is valid in these circumstances F is proportional to 


2/3 2/3 
R Gl Put hy e ho R 


Then, before peak intensity, de _ 


and after, dc _ -} Rc 
ater hoke(c 


Seen FOr Wo in theories 5b and 5c substitute Wo + W, where W is given by 


We 16W Vepiume 7? air» acting towards the centre of the array. Wy is 
the upward velocity of the plume ~ Thamas (4). 


Now drag force per unit height on plume is proportional to W_.Perimeter 4L . 
Assume this is equal, or at least proportional, to the buoyakicy force per 
unit height; and that this is a linear function of b' and the difference 
between the atmospheric lapse rate and the adiabatic lapse rate (where b! 

= total burning rate = sum (b) = R 

N.B. since buoyancy is proportional to L and drag to L, Wy is propor- 
tional to L and there must be an area effect. 


Buoyancy — depends 
on burning rote & 
lapse rate — |¥ Drag oc Wy. SL 


A180 Pyiume = CP air * (1-c)? £1 ame 


And , . Volume of excess air 
* " Li "+ volume.of flame 


With Volume of excess air = 4LnW - b!/ke 


and Volume of flame = b' , Yair . stoichiametric ratio 
ke Prlane 


Alternatively, assume the burning rate is a function of IW, + WI 


(cf theory 4a). In either case there is feedback on the burning rate. 


6 Thermal radiation feedback 
Burning rate b = % + by s 


Where S = solid angle of flame and burning fuel 
= S (piles) + S (other flames) 


Also S (flame) = S(b) - hence feedback. 
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Theories of Scaling 
7 Flame height scaling 


ja Flame height F proportional to linear dimensions L. 
Already falsified by trials - Craig Chandler 
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7b From Thoras (4), F/L is usually proportional to the cube root of 
reA. (note: this is a similar result to theory 9b below with n = 24). 
Other regimes are possible, but it is always a function of R@/L. 


Therefore F =a (w)2/3 (where a is a constant). 


Let piles be U cms square, with gaps of V ams between them. 


20a 
Then flame height for isolated pile = Fi =@ b it) 


And if the flames from N piles have merged, the new flame height, 
Fi «a v/3y'/3 (yg vy */3 


Now assume that isolated fires will join if F, >iV and that a nerged 


"damain" of N fires will persist. if Fo > 3 N1/2(U + V), where i, j are 
numerical constants, perhaps of the order of 1. (Thamas, Baldwin and 
Heselden's results (5) suggest that i is not constant, but is in the 
region of 5-7 in a Flambeau situation.) 


These correspond respectively to the ratio of the flame height to 
the gap, and to the dimensions of the damain. 


Assuming that these two conditions occur simultaneously, the ratio 
of fleme heights on either side of the discontinuity, (i.e. that with 
danains of 2 fires just not breaking down, divided by thet with isolated 
fires just not merging) will be:- 


2/2 4 (wey) 
ot a 


Thus three regimes are possible: isolated fires, fires merged to give 
shifting domsins of about N, or the whole array forming one giant flame. 


3) ONE GIANT FLAME 
In the partial merging regime, the fires will grow till 


Fe n'/2 (y 4 V) 5 is satisfied 


In this case b = wi/s (U + y)3/2 (4/a)3/2 


Thus N is proportional to Fe and to p* and Fo is proportional to pe. 


8 Fire intensity for a given fuel type depends only on fire loading and | 
is independent of area. 
x, High burning 


These proportionalities are stil) valid if F/L is any function of 


[p . rate 
it, 
LL, , 
fy J\ Lf Low burning Low burning 
rf I; \) rate rate 


Small area Small area large area 


High fire toad Low fire load Low fire lood 
9 Scaling of energy requirements. 


9a If all linear dimensions are to be doubled, then Power Needed 
increases by a factor of 8 and Burning rate/an.2 of ground (R) increases 
by 2. If Pressure difference does not vary then pressure gradient 
decreases by 2 and since distance air is to be accelerated increases by 
2 time needed to establish the airflow increases by 2. So minimum fire 
loading to establish the airflow increases by 4. 


9b As in theory 9a, but let power needed increase by oe 
Then K increases by 2n-2 
Time needed increases by 2 


Ninimaslcivebloadinusincreasagiby 2oe 


Tal 1: Ratios of quantities in 2 scaled bums 


Power needed L3 Ln 
Burning rate/am L | Ln-2 
Pressure Difference yo | io 
Pressure Gradient I! | L-1 
Time needed L L 
Minimum Fire Loading 12 Ln-1 


Theories of Weather and Fuel 


10 If initial temperature of wood «= rc, and breakdown temperature = 
T,°C, Burning rate is inversely proportional to (T, -T). 


7 
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11. If rain falls during a burn, heat production is reduced by 620 cals per 
gram of water falling on a pile. 


bo 
Heat 
Production | 
Rain to put fire out 
eras gmsisec 
——> Amount of rain 
12 If wood is wet, fire loading is reduced by 620 cals/gram of water. 


Effective fire 
loading 


Slope - 620 cals/gm 


<a Degree of wetness in wood 


13 If inert material is mixed with the fuel it has no effect on the heat 
generation; count only fuel for the fire losding. 


Adding inert 


F : material e.g.bricks — 
f/f 
7. \ no effect on the 
YS) ; ves { | burning rote 


Theories cf Shelter 


14 Theories of toxic gas penetration. 


CO Concentrotion C’ 


CO Concentration C 
Leak rate / air change: 


- et 


1ha Let C be the internal concentration of CO, and C' that outside. 
Then dC = p(Ct - C) 

dt 
and if D is the time integral of C over the whole period, and D' the 
integral of C', then by integrating the differential equation over the 
whole period, D = D' (equal dose theorem) 


Hence the shelter does no good, unless completely sealed (p< 0), 


Outside 


Inside 


Shaded areas are the same 


14b As in theory tha 
Put dC = pC! -q C - r (due to losses etc.) 
dt 
Then D = pD! - rt, 
q Tq 


where t, = total time with CO in the shelter 


2 


As a special case, CO will be lost by people breathing it. Assume 
that between than they consume a fraction f/unit time. 


therefore dC = p(C' -C) -£C 
dt 


thereforeq=p+f, r=0 


therefore D = pO! 
pef 


aye 


Some permanently lost 


14c As in theory 14a, but dose is not the criterion of harm; there is 
a threshold concentration Co. If this is reached at time t, in the shelter 
and ceases at time to: 


t 
Dae J. (C = 0) dt 
t 
2 
Pp J, (Grace 


Suppose that C' can be approximated by a constant value Q % for a time 
E/p. Then before this time, C » Q Cy (i -e 522) 


and after it C 22 Cy (63s 1) enh 


also t, = 1 log _Q_ 
P Q-1 


and t, = 1 log (2 (e#-1) 
Finally, D e Bt _ 110g Q __1 log {(Q - 1¢¢F-1)}) 
E Q@1° 


Note that D = O (perfect protection) when E < log _Q. 
Goa 


Outside 


Inside 


1kd Gravity stratification - Rasbash (6) 


Suppose gases enter the shelter through a window of area A with a 
speed W. The difference in density from the air in the shelter is Ap. 
Now if the window is at the bottom of the shelter there will be thorough 
turbulent mixing and theories 14a, b or c will apply; and also if the 
modified Richardson number _ Al/2 ap is less than 0.1. 

p 


RICHARDSON N2 <0-1 


Rapid Mixing 
—P by I6e 


Window area =A 


If it is more than 1, however, the gases in the shelter should be well 
stratified with the hot incoming gases lying on top. Assume the top part 
of the shelter fills up rapidly so that C = C'. The gases will then slowly 
diffuse downwards following the equation 


Se 1 2 
32° s “tt 


Where z is the vertical co-ordinate and s a diffusion coefficient. In the 
later stages, when C' starts dropping, assume that gravitational turnover 
will keep the maximum concentration down to C', 


RICHARDSON N2 >1 
Window area: A 


¥v 
Top part fills Slow downwerd Gravitational 
up rapidly diffusion turnover 
te=nl t:2 t:3 


lhe Assume there is "recovery" from CO poisoning so that the same dose 
over a longer time gives a smaller effect. If the recovery is linear 
this is equivalent to theory lic; if it is exponential (proportional to 
dose) or any other relationship, the protection can be evaluated 
numerically, knowing the dependence of C' on time ~ see theory 20. 


15 Theories of heat penetration. 


Etiee 


15a If there is a simple mixture of gases, then with ambient temperature T 


og 
corresponds to C in theory 14a. However, "Temperature Dose" is not the 
criterion of harm. 


, 


T 


Shaded areas equal 


To ; 
15b Allowing for the loss of heat to the shelter walls still leaves the 
mathematics of theory 14a valid. The constant p, however, is greatly 


reduced, by a factor heat capacity of air and walls. 
heat capacity of air | 


7 
|/ SHELTER 


AIR | 
=, EXTERIOR 
low heat . A 
Hightiheatz.copacltye Infinite heat capacity 


capacity ~ 


Heat Exchange 


15c Taking heat stroke as the limiting factor, the risk of this may be 
measured by the integral of T-T., where qT; is a critical temperature 


of about 32°C, The mathematics of theory 14c then follows. 


1g 
OUTSIDE 
Ts 
INSIDE 


t 


15d Skin burns occur when T is greater than T,, which is perhaps 400°C, 
If the probability of T being greater than T2 “at a specified moment is 
P, then divide the total exposure time t2 into effectively independent 
periods of duration to, The probability of a man getting burt at any 
time during t is then:- 


1 - exp ( - Pt,/t)) 


15e In theory 15d, to allow for the cambined effect of hot air and thermal 
radiation, let 1 cal/cm.2/sec. be equivalent to K°C, Tt can then be 
reduced by the appropriate amount. 


Threshold reduced 
* radiation 


0 to 2% 3f 4t, St. 6t) time 


15f Straight conduction of heat through the walls is probably not very 
important. An exact solution can be found by solving the differential 
equation of heat conduction in one dimension, given as a boundary con- 
dition the variation of the extemal temperature with time - see 
theory 20. 


16 Theories to calculate the leak rate. 


16a To calculate p in theory 14a, assume a room volume v,2 opposite open 
windows each of area A, and wind speed W. 


Then Py © AW/v 


Area A 


Windspeed W 
Volume v 


16b If there are pressure fluctuations of + B millibars, occurring with 
an average periodicity of t', 


Assume B is a function of the average windspeed only. 


16c Gravitational pumping - Heselden (7) 
If the height of the window fran top to bottom is H, then:~ 


Pp, =A (exe) 
3 4 


Low density in at top 


External density p - Ap 


Air 
density p 
High density out at 
bottom 
Volume v 


16d If there are wind, pressure and gravitational effects, 
Pi etBa te bo..* Ba 


NOTE that the wind, pressure and gravitational effects would all tend 
to be strongest when the fire is burning most fiercely. Hence p may 
tend to drop in the later stages of the fire, leading to retention of 
low concentrations of CO for longer than expected by theory tha. 


16e The graph of p against Wo (with no fire) is a characteristic of the 
structure depending only on geometry, and can be used to predict the 


value of p with a fire. 


No fire 
In fire 


Wo Wind Speed Sealed Crack Test Small Large 


Theories relating to measurements 
17. Gas Composition and Temperature. 


17a Air can be simply characterised by the proportion of burnt gases in it. 
Hence the temperature rise above ambient, the CO%, the CO02%, and the drop in 
02% are all proportional to one another, and a measurement of any one will 
enable all the others to be calculated. 


% of burnt gases 


° 
nn 


WA 
oer in the air 


17b Theory 17a is only true for a given type of burning; when this changes, 
the ratios are upset. 

18 Wood block instruments. 

18a The degree of darkening (due to scorching) of a piece of white wood 

or fibre board exposed in a fire area is a measure of the overall severity 
of fire conditions at that point and is directly related to the 

probability of ignition of cellulosic materials. 

18b The darkening is similarly related to the probability of skin burns. 


Wa&b 


Reflectivity Risk % 


Fire Severity 


19 If P is the probability of ignition of a flat model wooden structure, 
then the probability of ignition of a similar structure of M times the 
area is: 1-(1-P)M 


Probability of 
ignition 


Size of structure 


20 Fire effects vary with time in a standard fashion. 


20a Ali fire effects (thermal radiation, flame temperature, air temperature, 
soil tenperature, induced windspeed, rate of weight loss CO concentration, 
etc.) follow, on the average, identically the same time curve in all fires. 
This allows variations in a multiplying factor on the intensity scale, and 
also in a multiplying factor on the time scale (which will depend only on 
the burning rate, and so be the same for all effects in the same part of the 
seme fire), All temperatures are measured fran the initial temperature. 

The curves can be drawn on log:log paper. 


log Effect Rad" Burn} 


Rad” Burn? 


Temp Burnt 


Temp Burn 2 


ol 1 10 100 1000 Log. Time 


20b As in theory 20a, but each effect has a different characteristic 
shape of curve. 


1000, 


Log. Effect Red” Burnt 


100 J 
Rod" Burn2 


ea Paw 
Temp Burnt 
1 f oa Burn2 


“ ' 
o1 1 10 100 


1000 Log. Time 


21 Flame height (or length) F and angle from vertical 6 in a wind W with 
a turning rate b are given by:- 


F « F(b,W) 
6 = 6(b,W) 


These equations can be solved for b and W, given measurements of 
F and @ 


os 
Burning rote b 


2/3 : 
2la Fe a(uRy/? = a(b/u)~'~ as in theory 7b 
i.e.F is not affected by W or 0 


21b Measure K in mass units (gms lost/sec/horizontal sq. cm). Then gas 
must flow upwards with velocity R/p plame , and sideways with velocity W. 
Hence tan@ = W.Ppleme /R 


ee a - Speed = W 


Mass burning rote R gms env} sec 


Simple tests of some of the theories 


It should be fairly easy to investigate the variation of the burning 
rate or flame height with position in the array; some of the theories 
give quite different predictions on this:- 


2a Same everywhere. 


-~20- 


4a Maximum at or near edges falling steadily to minimum at centre - 
but may be modified by Wo 


me High on outer edge - all other piles low and same as each other. 
4c 


5a,5d,5e,5f All piles the same once well away from edges. 


~—— == 
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5b Progresses downwind through the array at constant speed - 
a wave motion. 


5c Maximum some distance in from the upwind edge. 


6 Low on outer edge - all other piles high and same as each other. 


oe ee Ee ee 
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